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Implementation  of  Wave  Dissipation 
by  Vegetation  in  STWAVE 

by  Mary  E.  Anderson  and  Jane  McKee  Smith 


PURPOSE:  This  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  describes  the 
implementation  of  wave  dissipation  by  vegetation  into  the  nearshore  spectral  wave  model 
STWAVE  (Massey  et  al.  2011;  Smith  2007;  Smith  et  al.  2001). 

INTRODUCTION:  The  influence  of  vegetation  on  coastal  hydrodynamics  is  a  relatively  new 
field,  with  a  body  of  literature  documenting  the  dissipation  of  wave  energy  by  coastal  vegetation 
developing  within  the  last  few  decades  (see  Anderson  et  al.  (2011)  for  a  summary). 
Unfortunately,  the  effect  of  vegetation  on  coastal  processes  and  hydrodynamics  is  not  fully 
implemented  in  many  numerical  models.  Standard  practice  in  nearshore  wave  propagation 
models,  including  STWAVE  and  SWAN,  is  to  account  for  energy  losses  due  to  vegetation  using 
bottom  friction  source  terms.  The  need  to  accurately  predict  coastal  hydrodynamics  in  the 
presence  of  natural  or  nature-based  features  has  led  to  an  increasing  demand  for  models  that 
better  capture  wave  interaction  with  vegetation.  Compounded  by  a  lack  of  technique  and 
guidance,  the  beneficial  effects  of  vegetation  are  often  neglected  in  the  analysis,  design,  and 
construction  of  coastal  protection. 

Theoretical  models  for  estimating  wave  dissipation  based  on  energy  conservation  were  initially 
proposed  for  monochromatic  waves  (Dalrymple  et  al.  1984),  with  later  expansions  to  narrow- 
banded  random  waves  (Mendez  and  Losada  2004).  One  noticeable  improvement  over  the  current 
bottom  friction  formulations  is  the  capability  to  describe  the  vegetation  itself.  The  declared 
vegetation  characteristics  are  the  following:  vegetation  height,  stem  diameter,  vegetation  density, 
and  a  bulk  drag  coefficient  calibrated  for  the  specific  plant  type  and  hydrodynamic  conditions. 
Calibration  of  this  bulk  drag  coefficient  accounts  for  many  processes  not  yet  fully  understood, 
such  as  plant  motion.  The  study  of  its  behavior  with  respect  to  wave  attenuation  is  currently 
ongoing.  The  random  wave  dissipation  model  proposed  by  Mendez  and  Losada  (2004)  is  the 
most  appropriate  for  inclusion  into  STWAVE  due  to  its  reasonable  representation  of  the  physical 
processes  and  feasibility  to  implement. 

FORMULATION:  Waves  propagating  through  vegetation  dissipate  energy  due  to  the  work 
carried  out  on  the  vegetation.  Assuming  the  validity  of  linear  wave  theory,  the  conservation  of 
energy  is  as  follows: 

dECv  P 

- -  =  — e 

dx 

where  E  is  wave  energy,  C,  is  group  velocity,  x  is  the  horizontal  distance  over  which  the  wave 
travels,  and  sv  is  the  time-averaged  rate  of  energy  dissipation  due  to  vegetation  per  unit  horizontal 
area.  Integrating  vertically  over  the  vegetation  height  (ls)  and  assuming  ev  is  only  a  function  of  the 
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horizontal  drag  force  (inertial  component  neglected),  the  definition  for  the  depth-integrated  and 
time-averaged  energy  dissipation  per  horizontal  area  for  a  vegetation  field  is  given  by 


e 


v 


Fxudz 


where  d  is  water  depth  and  z  is  the  vertical  dimension.  The  horizontal  force  per  unit  volume  (Fx) 
is  described  using  a  Morison-type  equation  (Morison  et  al.  1950): 


Fx=-pC„bvNu\u\ 


where  p  is  the  density  of  water,  C«  is  a  depth-averaged  drag  coefficient,  bv  is  stem  diameter,  N  is 
vegetation  density,  and  u  is  the  horizontal  fluid  velocity  due  to  wave  motion. 

Mendez  and  Losada  (2004)  modified  the  Dalrymple  et  al.  (1984)  fonnulation  by  using  a  Rayleigh 
distribution  to  describe  the  variation  in  wave  height,  yielding  the  following  for  random  waves: 


e 


v 


2*Jji 


PCD\N 


kg  f  sinh3  kls  +  3sinh  kls 
2(d)  3k  cosh3  kd 


H 


3 

rms 


where  g  is  the  acceleration  due  to  gravity,  k  is  wave  number,  co  is  wave  angular  frequency,  Ht7ns 
is  root-mean-square  wave  height,  and  CD  is  an  average  or  bulk  depth-averaged  drag  coefficient 
that  is  dependent  on  hydrodynamic  and  plant  characteristics. 

The  STWAVE  source  term  for  wave  damping  due  to  vegetation  (Sveg)  is  obtained  by  expanding 
the  Mendez  and  Losada  (2004)  formulation  for  sv  to  include  frequencies  and  directions  (Suzuki 
etal.  2011): 


s„,  =  -  jl g2cBbpN 

in  which  k  is  the  mean  wave  number,  co  is  the  mean  angular  frequency,  and  Etot  is  the  total 
wave  energy.  When  the  vegetation  height  exceeds  the  water  depth  (4  >  d),  the  vegetation  height 
is  dynamically  set  equal  to  the  water  depth  (4  =  d)  within  STWAVE. 

MODEL  INPUT.  The  primary  input  file  for  STWAVE  is  the  simulation  file  (*.sim),  and  it  is 
within  this  file  that  model  controls  are  defined  through  a  series  of  FORTRAN  namelists  that 
specify  parameters  and  run  options.  For  a  detailed  description,  see  the  STWAVE  v6.0  user’s 
manual  (Massey  et  al.  2011).  To  implement  vegetation  dissipation  into  STWAVE,  an  additional 
model  option  IVEG  was  required  in  the  std_parms  namelist.  A  description  of  the  IVEG  option  is 
provided  below. 

IVEG  =  Flag  to  exclude  vegetation  (IVEG  =  0,  default  option)  or  include  vegetation  (IVEG  = 
1  or  2).  Spatially  constant  vegetation  (IVEG  =  1)  requires  specification  of  vegetation 
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parameters  within  the  simulation  file  while  spatially  variable  vegetation  (IVEG  =  2) 
requires  an  input  file  that  specifies  the  vegetation  parameters  for  every  grid  cell.  The 
vegetation  parameters  required  for  both  options  are  the  following:  average  vegetation 
height,  average  stem  diameter,  number  of  plants  per  unit  horizontal  area,  and  bulk  drag 
coefficient. 

To  accommodate  the  spatially  constant  vegetation  option  (IVEG  =  1),  a  new  namelist  was  added 
called  const  veg.  The  vegetation  parameters  for  the  IVEG  =  1  option  must  be  defined  within  the 
const  veg  namelist  in  the  *.sim  file.  As  a  reminder,  all  FORTRAN  namelists  start  with  the 
ampersand  symbol  (&),  followed  by  the  namelist  name.  Variables  assigned  to  the  namelist  are 
then  listed  along  with  the  assigned  value.  The  end  of  the  namelist  is  indicated  by  the  slash 
symbol  (/).  The  const_veg  namelist  must  be  specified  before  the  optional  namelists,  which  are 
indicated  by  the  symbol  @  in  the  STWAVE  *.sim  file.  The  required  model  parameters  for  the 
new  const_veg  namelist  are  provided  in  Table  1. 


Table  1.  Model  parameters:  Spatially  constant  vegetation  -  const_veg  namelist. 

Parameter 

Type 

Definition 

veg_ls_const 

real  number 

#  =  average  vegetation  height  [m] 

veg_bv_const 

real  number 

#  =  average  stem  diameter  [m] 

veg_N_const 

real  number 

#  =  vegetation  density  [stems/rrk] 

veg Cd const 

real  number 

#  =  bulk  drag  coefficient  [-] 

For  the  spatially  variable  vegetation  option  (IVEG  =  2),  vegetation  parameters  for  each  grid  cell 
are  defined  using  an  external  input  file.  The  external  input  filename  is  specified  under  the 
input  file  namelist  in  the  *.sim  file  using  the  VEG  option.  The  values  must  be  provided  in 
column  format,  with  the  following  parameters  from  left  to  right:  average  vegetation  height, 
average  stem  diameter,  vegetation  density,  and  bulk  drag  coefficient.  The  format  is  the  same  as 
the  other  global  STWAVE  files  and  is  read  using  the  following  FORTRAN  algorithm: 

do  j  =  NJ,  1,-1 
do  i  =  1,  NI 

read  (10,  *)  veg_ls(ij),  veg_bv(i,j),  veg_N(i,j),  veg_Cd(i,j) 
enddo 
enddo 

Excerpts  of  the  required  changes  to  the  STWAVE  simulation  file  for  spatially  constant 
vegetation  (IVEG  =1)  and  spatially  variable  vegetation  (IVEG  =  2)  are  highlighted  in  Appendix 
A.  An  example  of  the  external  input  file  required  for  spatially  variable  vegetation  (IVEG  =  2)  is 
provided  in  Appendix  B. 

MODEL  VALIDATION 

Mendez  and  Losada  (2004).  Comparisons  of  STWAVE  to  the  analytical  dissipation  model  of 
Mendez  and  Losada  (2004)  were  completed  to  confirm  correct  implementation  of  the  source 
term.  The  random  wave  transformation  model  for  a  flat  bottom  was  proposed  by  Mendez  and 
Losada (2004) as 
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r,  1  „  ,  NR  ,  sinh3  k/s  +  3sinhk/s 
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where  HrmSi 0  is  the  root-mean-square  wave  height  incident  to  the  vegetation  (x  =  0).  The  root- 
mean-square  wave  height  is  related  to  the  significant  wave  height  ( Hmo )  by 


rms  V2 

STWAVE  simulations  were  carried  out  with  a  water  depth  (d)  of  1 .0  meter  (m),  peak  periods  (Tp) 
ranging  from  1.5  to  10.0  seconds  (s),  an  incident  significant  wave  height  (Hmo)  of  44.8  centimeters 
(cm)  at  the  input  boundary  (x  =  0),  and  a  spatially  constant  vegetation  field  measuring  100  m  in 
length.  The  STWAVE  domain  was  a  rectangular  grid  with  50  cells  in  the  alongshore  direction,  25 
cells  in  the  cross-shore  direction,  and  a  2.0  m  grid  resolution.  The  vegetation  characteristics  were 
varied  each  simulation  to  verify  the  model  accurately  captured  the  response  in  wave  height.  The 
model  parameters  for  all  simulations  are  given  in  Table  2.  Unidirectional,  narrow-banded  random 
waves  were  generated  using  Texel,  Marsen,  and  Arsole  (TMA)  shallow-water  spectra  with  a 
spectral  peak  enhancement  factor  (y)  of  10.0  and  a  directional  spreading  factor  (n)  (i.e.,  cosn6 )  of 
50.0.  Wave  breaking  and  bottom  friction  were  not  simulated  to  ensure  wave  energy  dissipation  was 
due  only  to  vegetation. 


Table  2.  Simulation  parameters  for  comparison  with  Mendez  and  Losada  (2004). 


Tp  [s] 

Is  [m] 

bv  [m] 

N  [stems/mz] 

CD  H  i 

1.5 

0.75 

200.0 

1.0 

3.0 

0.25 

0.02 

250.0 

0.2 

5.0 

1.25 

0.1 

10.0 

0.5 

8.0 

0.5 

0.004 

400.0 

0.35 

10.0 

1.0 

0.05 

100.0 

0.7 

Figure  1  demonstrates  the  comparisons  between  STWAVE  and  the  analytical  solution  of  Mendez 
and  Losada  (2004).  STWAVE  simulations  were  completed  both  in  half-plane  mode  (STWAVE- 
HP)  and  full-plane  mode  (STWAVE-FP)  to  ensure  the  solutions  were  consistent.  The  wave 
attenuation  trends  compare  well  with  each  other,  with  only  slight  differences  between  the  two 
STWAVE  modes  and  Mendez  and  Losada  (2004).  The  most  noticeable  difference  was  for  the 
shortest  peak  period  (Tp  =1.5  s).  These  inconsistencies  may  be  a  result  of  the  use  of  the  mean 
period  in  STWAVE  while  Mendez  and  Losada  (2004)  use  the  peak  period.  As  these  differences 
are  very  small,  the  wave  attenuation  proposed  by  Mendez  and  Losada  (2004)  for  unidirectional, 
nonbreaking  random  waves  can  be  closely  approximated  using  STWAVE. 
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Figure  1 .  Comparisons  of  wave  height  evolution  between  STWAVE  and  Mendez 
and  Losada  (2004). 
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SENSITVITY  TO  SPECTRAL  PARAMETERS 

Frequency.  The  sensitivity  of  the  STWAVE  solution  to  differences  in  the  peak  enhancement 
factor  (y)  is  investigated  by  comparing  TMA  spectra  o f  y  =  10.0,  considered  a  narrow  spectral 
value,  and  y  =  3.3,  a  common  default  value.  An  illustration  of  the  resulting  difference  in  spectral 
shape  is  provided  in  Figure  2.  A  larger  peak  enhancement  factor  (y)  results  in  a  narrower  and 
greater  concentration  of  energy  at  the  peak  frequency  compared  to  smaller  y  values. 


Figure  2.  TMA  spectra  with  y  =  10.0 
and  y  =  3.3  for  Tp  =  3.0  s. 


The  simulation  parameters  shown  in  Table  2  were  kept  constant  to  compare  results  to  those  of  the 
Mendez  and  Losada  (2004)  evaluation.  As  seen  in  Figure  3,  the  difference  in  wave  dissipation  due 
to  changes  in  peak  enhancement  factor  remained  negligible  throughout  the  range  of  considered 
periods. 

Direction.  Differences  in  wave  attenuation  due  to  changes  in  the  directional  spreading  factor  (n) 
are  also  investigated.  Input  spectra  for  both  simulations  was  TMA  spectra  with  y  =  3.3,  with  one 
spectrum  generated  with  n  =  50.0,  considered  a  narrow  spreading  value,  and  the  other  with  n  =  4.0, 
a  standard  value.  Again,  all  other  simulation  parameters  remained  constant.  Energy  is  distributed 
over  a  narrower  angle  band  for  larger  values  of  n,  as  can  be  seen  in  Figure  4. 

The  difference  in  wave  dissipation  for  a  spatially  constant  vegetation  field  is  shown  in  Figure  5. 
Little  difference  in  wave  evolution  due  to  changes  in  directional  spreading  was  observed  along 
the  center-axis,  with  these  differences  remaining  small  throughout  the  range  of  tested  periods, 
although  the  n  =  4.0  solution  is  consistently  slightly  lower  than  n  =  50.0. 

The  difference  in  wave  dissipation  considering  a  patch  of  vegetation  smaller  than  the 
computational  grid  was  also  investigated.  These  computations  were  completed  on  a  500  m  long 
by  300  m  wide  rectangular  grid,  with  a  smaller  vegetation  patch  along  the  center-axis.  The 
vegetation  patch  was  204  m  long,  120  m  wide,  and  located  approximately  134  m  from  the 
offshore  boundary.  The  water  depth  remained  constant  at  d  =  1.0  m  throughout  the  domain.  The 
modeled  incident  significant  wave  height  and  peak  period  were  Hmo  =  44.8  cm  and  Tp  =  5.0  s. 
Local  vegetation  was  constant  at  ls  =  60.0  cm,  bv  =  5.0  millimeters  (mm),  N  =  250.0  stems/m2, 
and  Cd  =  1.0.  Simulations  were  completed  using  STWAVE -HP. 
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Figure  3.  Wave  height  evolution  for  TMA  spectra  with  y  =  10.0 
and  y  =  3.3  and  n  =  50.0. 


Figure  4.  TMA  spectra  with  n  =  50.0  and  n  =  4.0  for  Tp  =  5.0  s. 


Figure  6  and  Figure  7  show  that  the  most  noticeable  difference  in  wave  height  is  located 
downstream  of  the  vegetation  patch.  Reduced  wave  height  is  seen  behind  the  patch  for  n  =  50.0  as 
the  waves  are  impeded  and  limited  in  direction.  However,  waves  are  able  to  propagate  and  refonn 
downstream  of  the  vegetation  for  n  =  4.0  as  wave  energy  is  distributed  amongst  a  broader  range  of 
directions.  This  wave  behavior  is  similar  to  that  observed  for  irregular  wave  interactions  with 
structures. 
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Figure  5.  Wave  height  evolution  for  TMA  spectra  with  y  =  3.3  and 
n  =  50.0  and  n  =  4.0. 
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>% 


x  [m] 


Figure  6.  Wave  interaction  with  vegetation  patch  for  TMA  spectra  with  y  = 
3.3  and  n  =  50.0  (top)  and  n  =  4.0  (bottom). 
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Figure  7.  Wave  height  evolution  along 

center-axis  for  TMA  spectra  with  y 
=  3.3  and  n  =  50.0  and  n  =  4.0. 

APPLICATION  TO  LABORATORY  EXPERIMENTS:  Two  bottom  friction  dissipation 
fonnulations  are  available  in  STWAVE,  JONSWAP  and  Manning’s  n.  The  vegetation  source 
term  was  combined  directly  with  the  bottom  friction  terms  in  order  for  dissipation  to  be  a 
function  of  the  individual  processes  (bottom  friction  or  vegetation)  or  both  processes 
concurrently  (bottom  friction  and  vegetation).  To  test  this  formulation,  the  STWAVE  model  with 
vegetation  was  applied  to  laboratory  measurements  of  wave  attenuation  conducted  by  Anderson 
and  Smith  (2014).  The  experiment  was  carried  out  in  a  63.4  m  long,  1.5  m  wide,  and  1.5  m  deep 
wave  flume.  The  vegetation  zone  measured  9.8  m  long  and  was  populated  with  idealized 
vegetation  (Figure  8).  The  idealized  vegetation  was  constructed  from  6.4  mm  diameter  flexible 
tubing  measuring  41.5  cm  long.  Two  densities,  200  and  400  stems/m2,  were  tested  for  21 
irregular  wave  conditions  with  varied  wave  height,  peak  period,  and  water  depth.  Input  spectra 
for  STWAVE  were  generated  from  the  measured  spectra,  with  all  energy  focused  in  the  0-degree 
angle  band  as  waves  were  unidirectional. 


Figure  8.  Model  setup  for  Anderson  and  Smith  (2014). 


To  test  the  implementation  of  the  dissipation  fonnulation  with  the  measured  data,  wave 
dissipation  was  considered  in  the  following  form: 
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where  the  first  term  on  the  right-hand  side  of  the  equation,  shown  in  bold,  is  the  Manning’s  n 
bottom  friction  source  tenn,  and  the  second  term  is  the  vegetation  source  tenn.  The  background 
dissipation  due  to  friction  in  the  flume  (Cbj),  was  first  solved  for  by  assuming  no  vegetation: 


Cbf=-  In 


K' 

g 

sinh  kh' 

A 

{H2o\ 

^ rms 

.  co  , 

Ax 

where  the  root-mean-square  horizontal  wave  velocity  (urms),  was  defined  at  the  bottom  and  was 
calculated  as 


H0  co 
2 V2  sinh  kh 


Cbf  was  calculated  using  the  unimpeded  control  runs  (no  vegetation),  with  Hmo  measured  at  the 
gauge  immediately  incident  to  the  vegetation  test  section  and  at  the  last  gauge  corresponding  to 
Ho  and  H2,  respectively.  The  location  of  these  gauges  of  interest  is  identified  in  Figure  8. 

After  estimating  the  background  friction,  the  drag  coefficient  due  to  vegetation  was  estimated: 
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Again,  Ho  is  measured  at  the  gauge  immediately  upstream  of  the  vegetation  test  section,  and  1 12 
is  measured  at  the  last  gauge  (downstream  end).  The  input  into  STWAVE,  Manning’s  n,  was 
obtained  by  rearranging  Cbf. 
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bf  ~ 


h1/3 

bf'1 
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Figure  9  compares  the  nondimensional  wave  evolution  of  STWAVE  with  a  subsample  of  the 
flume  experiments.  The  top-left  panel  shows  the  wave  attenuation  with  a  water  depth  of  d  =  30.5 
cm,  an  incident  significant  wave  height  of  Ho  =  1 1.3  cm,  a  peak  wave  period  of  T„  =1.25  s,  and  a 
stem  density  of  N  =  400.0  stems/m'.  The  top-right  panel  shows  the  wave  attenuation  with  a  water 
depth  of  d  =  45.7  cm,  an  incident  significant  wave  height  of  H0  =  19.2  cm,  a  peak  wave  period  of 
Tp  =  2.0  s,  and  a  stem  density  of  N  =  200.0  stems/nr.  The  bottom  panel  shows  the  wave 
attenuation  with  a  water  depth  of  d=  53.3  cm,  an  incident  significant  wave  height  of  Ho  =  11.1  cm, 
a  peak  wave  period  of  Tp=  1.5  s,  and  a  stem  density  of  A  =  200.0  stems/m".  Overall,  STWAVE  is 
able  to  replicate  the  wave  evolution  trend  of  the  experiments  well,  with  a  goodness  of  fit 
coefficient  ( R~ )  exceeding  0.90  for  all  comparisons.  A  major  assumption  in  implementation  of 
vegetation  is  the  linear  summation  of  bottom  friction  and  vegetation  losses.  This  model  validation 
supports  that  this  approximation  is  appropriate. 
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Figure  9.  Wave  height  comparisons  of  STWAVE  and  the  laboratory  results  of 

Anderson  and  Smith  (2014):  (top  left)  d  =  30.5  cm,  H0=  1 1.3  cm,  Tp  =  1 .25 
s,  N  =  400.0  stems/m  ,  (top  right)  d  =  45.7  cm,  H0  =  19.2  cm,  Tp  =  2.0  s,  N 
=  200.0  stems/m2;  (bottom)  d  =  53.3  cm,  H0=  11.1  cm,  Tp  =  1.5  s,  N  = 

200.0  stems/m2. 

CONCLUSIONS:  The  analysis,  design,  and  construction  of  coastal  protection  often  neglect  the 
beneficial  effects  of  natural  or  nature-based  features  because  insufficient  methods  are  available 
to  capture  those  benefits.  In  this  technical  note,  random  wave  dissipation  by  vegetation  is 
implemented  into  the  phase-averaged  nearshore  wave  model  STWAVE.  The  vegetation 
dissipation  source  term  is  first  validated  by  comparing  model  results  to  the  solution  of  Mendez 
and  Losada  (2004)  for  unidirectional,  nonbreaking  waves.  Following  validation,  the  sensitivity  of 
the  model  solution  to  energy  distribution  in  frequency  and  direction  is  investigated.  The 
difference  in  wave  height  as  a  result  of  frequency  spread  was  minute  for  the  range  of  considered 
periods.  The  wave  height  evolution  along  the  center  line  of  a  domain  with  spatially  constant 
vegetation  was  also  insensitive  to  differences  in  directional  spreading.  However,  the  global  wave 
solution  changes  considerably  with  directional  spreading  when  considering  a  local  vegetation 
patch,  particularly  downstream  of  the  vegetation.  Finally,  combining  the  vegetation  source  term 
with  the  already  existing  bottom  friction  terms  in  STWAVE  was  tested  using  the  flume 
experiments  of  Anderson  and  Smith  (2014).  By  first  solving  for  the  bottom  friction  coefficient 
and  then  the  bulk  drag  coefficient,  a  good  prediction  of  wave  dissipation  was  obtained.  Since  the 
bulk  drag  coefficient  accounts  for  many  assumptions  and  processes  not  yet  fully  understood, 
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calibration  of  the  bulk  drag  coefficient  is  essential  to  obtain  accurate  results,  and  its  dependence 
on  hydrodynamics  and  plant  biomechanics  requires  additional  investigation. 

SYMBOLS 

k  mean  wave  number 

co  mean  wave  angular  frequency 

bv  average  stem  diameter 

Cbf  bottom  friction  coefficient 

CD  average  or  bulk  depth-averaged  drag  coefficient 

Cd  depth-averaged  drag  coefficient 

Cg  group  speed 

d  water  depth 

E  wave  energy 

Etot  total  wave  energy 

Fx  horizontal  force  per  unit  volume  on  a  stem  array 

g  acceleration  due  to  gravity 

H  local  significant  wave  height 

Hmo  significant  wave  height 

EI^  root-mean-square  wave  height 

k  wave  number 

4  average  vegetation  height 

n  directional  spreading  factor 

N  vegetation  density 

tibf  Manning’s  n 

Sveg  wave  damping  due  to  vegetation  source  term 

Tp  peak  period 

u  horizontal  wave  particle  velocity 

Urms  root-mean-square  horizontal  wave  particle  velocity 

x  horizontal  distance 

y  spectral  peak  enhancement  factor 

sv  time-averaged  rate  of  energy  dissipation  per  unit  horizontal  area  due  to  vegetation 

p  water  density 

co  wave  angular  frequency 

z  vertical  dimension 

ADDITIONAL  INFORMATION:  This  CHETN  was  prepared  as  part  of  the  Wave  Dissipation  by 
Vegetation  for  Coastal  Protection  work  unit  in  the  Flood  and  Coastal  Systems  R&D  Program  and 
was  written  by  Mary  E.  Anderson  (Mary. Anderson®, usace.army.mil)  and  Jane  Smith 
(Jane.M.  Smith®, us  ace.  army,  mil )  of  the  U.S.  Anny  Engineer  Research  and  Development  Center 
(ERDC),  Coastal  and  Hydraulics  Laboratory  (CHL).  The  Program  Manager  is  Dr.  Cary  Talbot, 
and  the  Technical  Director  is  William  Curtis.  This  CHETN  should  be  cited  as  follows: 

Anderson,  M.  E.,  and  J.  M.  Smith.  2015.  Implementation  of  wave  dissipation  by 
vegetation  in  STWAVE.  ERDC/CHL  CHETN-I-85.  Vicksburg,  MS:  U.S.  Army 
Engineer  Research  and  Development  Center,  http://chl.erdc.usace.army.mil/chetn. 
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APPENDIX  A 

Example  STWAVE  simulation  file  with  IVEG  =  1  (constant  vegetation) 

#  STWAVESIMFILE 

#  written  from  SMS  11.1.5  64-bit 

# 

# 

mmmmmmmmmmmmmmmmmmmmmmm 

# 

#  Standard  Input  Section 

# 

&std_parms 
iplane  =  0, 
iprp  =  1, 
icur  =  0, 
ibreak  =  0, 
irs  =  1 , 
nselct  =  50, 
nnest  =  0, 
nstations  =  0, 
ibnd  =  0, 
ifric  =  0, 
idepopt  =  1 , 
isurge  =  0, 
iwind  =  0, 
i_bcl  =  2, 
i_bc2  =  3, 
i_bc3  =  0, 
i_bc4  =  3, 
iveg  =  1 
/ 

# 

#  Constant  Vegetation 

# 

&const_veg 
veglsconst  =  1.5 
vegbvconst  =  0.04, 
vegNconst  =  200, 
vegCdconst  =  0.5 
/ 
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Example  STWAVE  simulation  file  with  IVEG  =  2  (VARIABLE  vegetation) 

#  STWAVESIMFILE 

#  written  from  SMS  11.1.5  64-bit 

# 

# 

############################################## 

# 

#  Standard  Input  Section 

# 

&std_parms 
iplane  =  0, 
iprp  =  1, 
icur  =  0, 
ibreak  =  0, 
irs  =  1 , 
nselct  =  50, 
nnest  =  0, 
nstations  =  0, 
ibnd  =  0, 
ifric  =  0, 
idepopt  =  1 , 
isurge  =  0, 
iwind  =  0, 
i_bcl  =  2, 
i_bc2  =  3, 
i_bc3  =  0, 
i_bc4  =  3, 
iveg  =  2 
/ 

# 

#  Input  Files  Section 

# 

&input_fdes 

DEP  =  "vegtest.dep", 

SPEC  =  "vegtest.eng", 

VEG  =  "veg.in" 

/ 
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APPENDIX  B 

EXAMPLE  STWAVE  GLOBAL  VEG  INPUT  FILE  FOR  IVEG  =  2  (VARIABLE  VEGETATION) 

#  STWAVESPATIALDATASET 

# 

&DataDims 
DataType  =  0, 

NumRecs  =  1 , 

NumFlds  =  4, 

NI  =  50, 

NJ  =  15, 

DX  =  2.0, 

DY  =  2.0, 

GridName  =  "vegtest.sim" 

/ 

# 

&Dataset 

FldName(l)  =  "Average  vegetation  height", 

FldUnits(l)  =  "m", 

FldName(2)  =  "Average  vegetation  diameter", 

FldUnits(2)  =  "m", 

FldName(3)  =  "Vegetation  density", 

FldUnits(3)  =  "stems/mA2", 

FldName(4)  =  "Bulk  drag  coefficient", 

FldUnits(4)  = 

Reclnc  =  1 

/ 

IDD  veg 


1.5 

0.04 

200 

0.5 

1.5 

0.04 

200 

0.5 

1.5 

0.04 

200 

0.5 

1.5 

0.04 

200 

0.5 

1.5 

0.04 

200 

0.5 

NOTE:  The  contents  of  this  technical  note  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an  official 
endorsement  or  approval  of  the  use  of  such  products. 
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